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ABSTRACT 

The first generation of stars is thought to have formed in low-mass halos with T v ; r < 10 4 K where H2 
cooling is paramount. However, the efficiency of H2 formation and cooling in these halos may have been 
severely limited by feedback processes. In this paper we investigate the radiative cooling and collapse 
of halos with virial temperatures T V i r > 10 4 K, i.e. those that can cool in the absence of H2 via neutral 
atomic lines. The evolution of these halos differs from their less massive counterparts. Efficient atomic 
line radiation allows rapid cooling to ~ 8000 K; subsequently the gas can contract nearly isothermally 
at this temperature. In the absence of H2 molecules, the gas would likely settle into a locally stable 
disk and only disks with unusually low spin would be unstable. However, we find that the initial atomic 
line cooling leaves a large, out-of-equilibrium residual free electron fraction. This allows the molecular 
fraction to build up to a universal value of xh 2 ~ 10 -3 , almost independently of initial density and 
temperature. We show that this is a non-equilibrium freeze-out value that can be understood in terms 
of timescale arguments. Unlike in less massive halos, H2 formation and cooling is largely impervious 
to feedback from external UV fields, due to the high initial densities achieved by atomic cooling. The 
newly formed molecules cool the gas further to ~ 100 K, and allow the gas to fragment on scales of a 
few xlOO Mq. We investigate the importance of various feedback effects such as H2 photodissociation 
from internal UV fields and radiation pressure due to Lya photon trapping, which are likely to regulate 
the efficiency of star formation. 

Subject headings: cosmology: theory - early universe - galaxies: formation - molecular processes 



1. INTRODUCTION 

The first generation of stars and/or quasars must have 
formed out of gas with primordial composition dictated 
by Big Bang nucleosynthesis. Only after the first genera- 
tion of stars explode as supernovae and achieve widespread 
metal pollution is metal line cooling possible. Thus, the 
star formation efficiency and the initial mass function 
(IMF) of Population III stars is likely to be set by the 
physical regime in which metal-free cooling can take place. 
In recent years, the cooling of metal-free gas in the first 
halos that are able to cool within a Hubble time i coo i < tn, 
has come under intensive study (see Abel & Haiman 2000 
for a recent review). Such halos have low virial tempera- 
tures, T v ir < 10 4 K; at these temperatures H 2 is formed 
via gas-phase processes such as H + c~^H~+7, fol- 
lowed by H~ + H — > H2 + e~. Detailed numerical simula- 
tions have shown convergence toward a regime T ~ 200 K , 
n ~ 10 4 cm -3 , dictated by the thermodynamic properties 
of H 2 (Abel et al. 2000; Bromm et al. 2001a), which allows 
gas fragmentation into clumps of mass 10 2 — 10 3 M©. 

However, as the first stars begin to shine, they emit 
photo-dissociating radiation in the Lyman- Werner bands 
(LW; 11.2-13.6 eV) to which the universe is optically thin, 
and further H2 formation and cooling can be suppressed 
both by external (Haiman, Rees & Loeb 1997; Haiman, 
Abel & Rees 2000; Ciardi et al 2000; Machacek et al. 



2001; Ricotti, Gnedin & Shull 2001a, 2001b) and inter- 
nal (Omukai k, Nishi 1999; Glover & Brand 2000) sources 
of UV radiation. This may be partially alleviated if X-rays 
with energies <; 1 keV (whose mean free path in a uni- 
form IGM exceeds the Hubble length) can boost the free 
electron fraction, and thus the H2 formation rate in cool- 
ing clumps (Haiman, Rees & Loeb 1996; Haiman, Abel & 
Rees 2000; Oh 2001); note that such effects are unlikely to 
be important in the low density environment of the IGM 
(Venkatesan, Giroux, & Shull, 2001). Alternatively, pro- 
togalactic shocks can produce substantial H 2 (Shapiro & 
Kang 1987; Ferrara 1998), or H2 may be generated in front 
of HII regions or in relic HII regions (Ricotti, Gnedin & 
Shull 2001a, 2001b), or the optical depth in intergalactic 
H2 may reduce the strength of the external dissociating 
radiation field (Ricotti, Gnedin & Shull 2000). However, 
while these effects may promote H2 formation in the pres- 
ence of an external UV background radiation field, they are 
unlikely to balance the strong internal feedback provided 
by stars within a galaxy: in the worst case scenario, only 
one star per halo can be formed (e.g., Madau & Rees 2001). 
It is therefore often believed that efficient and widespread 
star (and/or quasar black hole) formation capable of reion- 
izing the universe must await the collapse of halos with 
T vir > 10 4 K, or M ha i > 10 8 [(1 + z)/ll]- 3 / 2 M (Ostriker 
& Gnedin 1996, Haiman, Rees & Loeb 1997, Ciardi et 
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al 2000, Abel & Haiman 2000, Ricotti, Gnedin & Shull 
2001a, 2001b). These halos do not rely on the presence of 
H2 molecules, since they can cool via recombination and 
collisional excitation of neutral atoms. 

The goal of this paper is to critically examine the pre- 
vailing assumptions of efficient gas cooling and star forma- 
tion in metal free halos with T V - 1T > 10 4 K. To date, the 
details of gas cooling and chemistry in these halos have not 
been studied with the same care and attention devoted 
to lower mass halos. Since the majority of stars and/or 
quasars which reionized the universe, and polluted the in- 
tergalactic medium (IGM) with metals, are expected to 
form in such halos, it is important to study the evolution 
of the gas in these halos in detail. 

It is generally taken for granted that Lya cooling of 
neutral atomic hydrogen allows rapid contraction of gas 
until it becomes self-gravitating at the center of the po- 
tential well and fragments to form stars (ideas tracing 
back to Rees & Ostriker 1977; White & Rees 1978). How- 
ever, it is also known that significant contraction is re- 
quired for this fragmentation to result in stellar-mass 
fragments. If only Lya cooling operates, the gas re- 
mains at a temperature of ~ 10 4 K due to the sharp cut- 
off in the (equilibrium) cooling function and the Jeans 
mass is exceedingly high, even at high densities: Mj ss 
10 8 (T/10 4 K) 3 / 2 (n/l cm- 3 )- 1 / 2 M . Even if the gas be- 
came self-gravitating, unless the gas can contract to ex- 
tremely high densities, n > 10 12 cm -3 , fragmentation to 
lower Jeans masses M ~ 100M Q cannot proceed. How- 
ever, the gas cannot cool to arbitrarily high densities, but 
eventually must form a rotationally supported disk; at the 
temperatures allowed by Lya cooling, we shall show that 
the majority of such disks are locally gravitationally sta- 
ble. 

Thus, an additional coolant is needed to lower the gas 
temperature, both to ensure gravitational instability and 
to lower the Jeans mass by several orders of magnitude. In 
the absence of metals, H2 formation and cooling is there- 
fore still critical to star formation in T V - 1T > 10 4 K halos, 
and cannot be ignored. In this paper, we study H2 for- 
mation in shock-heated gas, and show that a universal H2 
fraction of ih 2 ~ 10~ 3 forms in gas that cools from an 
initial temperature of T > 10 4 K. Similar behavior has al- 
ready been noted in previous studies of pregalactic shocks 
(Shapiro & Kang 1987, Kang & Shapiro 1992). Here we 
study the non-equilibrium H2 formation in the T v ; r > 10 4 K 
halos of interest, examine the robustness of this mechanism 
to variations in density, temperature, and radiation field, 
and show that the asymptotic abundance can be under- 
stood in terms of timescale arguments. These arguments 
reveal that over a wide range of densities, the cooling gas 
follows a universal track in the (x e ,T) plane. 

This paper is organized as follows. In §|[ we study 
the equilibrium structure of isothermal disks embedded in 
dark matter halos, and show that H2 cooling is needed to 
promote gravitational instability in most disks. In §|[ we 
use semi-analytic methods and non-equilibrium chemistry 
to investigate H2 formation and radiative cooling in halos 
with virial temperatures T v - n > 10 K, and argue that the 
H2 abundance builds up to a universal value of xn 2 ~ 10 -3 
under most realistic conditions. In §|, we study the ef- 
fects of H2 destruction by internal and external sources of 
UV radiation and show that feedback processes are much 



less efficient in T V1I > 10 K halos than in their smaller 
counterparts, primarily because in the larger halos, gas 
can be compressed to high densities by initial atomic line 
cooling. We also examine whether opacity and radiation 
pressure effects can halt the collapse or fragmentation. Al- 
though we conclude that this is unlikely, we argue that it 
could affect the efficiency of star formation. In we 
summarize our conclusions and the implications of this 
work. In all numerical estimates, we assume a cold dark 
matter cosmology with a cosmological constant (ACDM) 
(n„,n A> n b ft a ,/i > £r 8A - 1 ) = (0.3,0.7,0.019,0.7,1.0) (see, 
e.g., Bahcall et al. 1999 for a review of these choices). 

2. DISK FORMATION AND GRAVITATIONAL INSTABILITY 

Let us first consider the cooling and collapse of an ini- 
tially spherical configuration of gas in a typical halo with 
T v ir > 10 4 K. The discussion presented in this section 
serves two purposes: it will highlight the importance of 
H2 molecules for the larger halos, and will also yield the 
physically appropriate range of density, temperature and 
ionization fractions under which to consider H2 formation 
and cooling in later sections. 

The well-known condition for runaway contraction of 
gas is t coo \ < idyn < £ sc (see Rees & Ostriker 1977), 
which is typically satisfied for hydrogen atomic line cool- 
ing. Since cooling due to Lya has a very sharp cutoff at 
T < 10 4 K, the cooling time is a very sensitive function of 
temperature, increasing rapidly by several orders of mag- 
nitude in a very narrow temperature range below ~ 10 4 K 
(Spitzer 1978). This is because the free electron fraction 
able to excite atomic line cooling drops very rapidly in 
this temperature range. Thus, the condition t coo \ < idyn 
implies that the collapse is nearly isothermal. If the gas 
cools below ~ 10 4 K, the gas recombines, the cooling time 
rapidly increases until i coo i > idyn, and further contrac- 
tion of the gas is close to adiabatic, with the gas heat- 
ing up due to the contraction until the hydrogen atoms 
are collisionally re-ionized and the condition t coo i < idyn 
is satisfied once again. In the absence of any other ef- 
fects, this thermostatic mechanism would allow the gas in 
a T v ; r > 10 4 K halo to cool and contract to arbitrarily high 
densities. In practice, however, the gas has some initial an- 
gular momentum, and must therefore eventually become 
rotationally supported. After a contraction by a factor 
of A" 1 - 20 in radius (where A = J\E\^ 2 /GM 5 / 2 is the 
spin-parameter, and J, E, and M are the total angular 
momentum, energy and mass of the halo), this results in 
rotationally supported disk at the center of the halo (Mo, 
Mao & White 1998, hereafter MMW). It is possible that 
some fragmentation takes place as the gas collapses toward 
a disk, however, this is likely to be inefficient due to the 
slow growth of density fluctuations in a rapidly contract- 
ing medium (e.g., Kashlinsky & Rees 1983). We therefore 
assume that most fragmentation must take place in the 
disk itself. We shall show that if only atomic line cooling 
operates, the disk will be stable to fragmentation in the 
majority of cases. For fragmentation to proceed, an ad- 
ditional cooling channel (such as H2) is required to allow 
cooling below ~ 10 4 K. 

For simplicity, we shall follow MMW, and assume that 
the gas settles to an isothermal, exponential disk, with 
gas temperature Tgas, embedded in a halo of virial tem- 
perature T v ir with a Navarro, Frenk & White (1997, here- 
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after NFW) dark matter density profile. We begin by 
listing the characteristic properties of such disks, which 
will be relevant for our later studies of H2 formation 
and self shielding. Let us assume that baryons make 
up the universal mass fraction f^b/^m of the halo, of 
which some fraction /<j have collapsed into the disk, i.e. 
-/Wdi sk = TO d Mhai = /d(^b/^m)Mh a i . We also assume 
that the angular momentum of the disk J<j is some frac- 
tion jd of the halo angular momentum J, i.e. Jd = jdJ- 
Henceforth, we shall assume that the specific angular mo- 
mentum of the disk is similar to that of the halo, and thus 
m d = j d - 

The assumption that the baryons preserve their specific 
angular momentum during collapse results in a good fit to 
the observed size distribution of galactic disks (Mo, Mao 
& White 1998). Detailed numerical simulations (Navarro 
& White 1993, Navarro & Steinmetz 2000) have not sup- 
ported this simple model, but instead produced signifi- 
cantly smaller disks, due to the transfer of angular mo- 
mentum from the gas to the dark matter during the highly 
inhomogeneous collapse. These simulations, however, also 
fail to produce the sizes and properties of observed galaxy- 
sized disks. Inclusion of suppression of cooling until late 
times (Weil, Eke & Efstathiou 1998) or supernovae feed- 
back (Thacker & Couchman 2001) reduces the discrep- 
ancy, but the issue has yet to be conclusively resolved. 
We therefore simply extend successful semi-analytic mod- 
els of disk formation at low redshift to higher rcdshift, and 
note that future numerical simulations of such halos may 
not in fact produce such disks. 

The hydrogen number density at radius r and at vertical 
height z in an isothermal exponential disk of radial scale 
length i?d is given by (Spitzer 1942) 



n(r, z) — n Q exp 
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where n is the central density, z is the vertical scale 
height of the disk at radius r, 

z = ^ (o) 

(4TrGtim H n e- 2r / R *y/ 2 1 ' 

c s is the sound speed of the gas, and /1 = 0.6 is the mean 
molecular weight (we choose the definition of the disk 
scale length Rd to conform to the customary assumption 
that the surface density has an exponential profile: S oc 
exp(— r/Rd))- For a disk in a halo with spin parameter A, 
if we assume the baryons conserve their specific angular 
momentum when they collapse, then the disk scale length 

r 200j 
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is given by R d = 2- 1 ^(j d /m d )\r 2 oa.fc 1 .fE 
where r 2 oo ~ ?Vir is the radius that encloses a mean inte- 
rior mass density of 200p cr i t , / c (c) and /r(A, c, m d ,j d ) are 
dimcnsionless functions of order unity (MMW), and c is 
the dimcnsionless concentration parameter (NFW) . 

The central number density of the gas is obtained by 
setting J dz J 2irr drfj,m p n(r, z) = Mdi s k, which yields: 
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When considering the characteristic densities for H2 for- 
mation, it will be sufficient to consider densities at most 
an order of magnitude below the central density: gas with 
n > 0.1n o comprises > 50% of the mass of the disk. Be- 
cause hydrostatic support is only relevant in the vertical 

direction, as the gas cools the disk becomes thinner, with 

1/2 

a reduced scale height z a cx T g i s . The vertical column 
density of gas as a function of radius is: 



N m (r) = \/2n z exp(-2r/i?d) 
= 3 x 10 23 exp(-r/i? d ) 
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which is sufficient for self-shielding of the gas against both 
ionizing UV radiation and H2 dissociating radiation in the 
LW bands to become important. Note that the column 
density is independent of the gas temperature. 

We now consider the conditions for gravitational insta- 
bility of the disk. In computing the rotation curve V(r), 
we use the formalism of MMW, which takes into account 
the contraction induced in the inner regions of the halo 
by the cooling and formation of the disk. This is done 
by assuming the disk is assembled slowly and the angular 
momentum of dark matter particles is an adiabatic invari- 
ant (Blumenthal et al. 1986, Flores et al. 1993). For 
the disk to be locally gravitationally unstable despite the 
stabilizing effects of tidal shears and pressure forces, we 
require the Toomre parameter Q < 1, where (e.g. Binncy 
& Tremaine 1987) 

(5) 
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£ is the disk surface mass density, and k = lAl(V/r)(l + 
dlnVydlnr) 1 / 2 is the epicyclic frequency. Regions in local 
disk galaxies where Q > 1 are observationally associated 
with very little star formation, indicating that the Toomre 
criterion is obeyed remarkably well (Kennicutt 1989). For 
our purposes, if Q > 1 everywhere throughout the disk, 
it is gravitationally stable and we assume no star forma- 
tion takes place. Disks with high spin parameters have low 
surface densities and satisfy this criterion. For any given 
disk-halo system, we can calculate a critical spin param- 
eter A cr it for which the disk is marginally stable. Here 
we define A cr ;t by the requirement that Q attains a mini- 
mum value of Q = 1 at least at one position in the disk. 
The fraction of halos that remain dark can therefore be 
found by integrating over the spin parameter distribution, 
/dark = J~ rlt p(A)dA, where p(X) is given by 



p(X)d\ = 



*A(27r)V: 



;CXp 



ln 2 (A/A) A dX 
J A 



2a 2 



(6) 



^0.057 v 10 

2 In another context, it was proposed that similarly high spin dark halos 
ct al. 1997). 



with A = 0.05 and o\ = 0.5 (e.g. Barnes & Efstathiou 
1987, Warren et al. 1992). 2 Our use of the Toomre cri- 
terion to characterize gravitational instability depends on 

might be detected from their gravitational lensing signature (Jimenez 
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the assumption that the disk is thin, z a <C Rd] in this 
case gravitational instability is insensitive to the vertical 
structure of the disk (Goldreich & Lynden-Bell 1965). The 
ratio of scale heights is given by 

* = 1.5x10- (JSL-) ( T ™ V 1 (A) 
R d V8000K7 \5 x 10 4 Ky \0.05 J \0.5 J 

,(7) 

and thus z D <C i?d is generally satisfied, except in low virial 
temperature halos where T v ; r ~ T gas . 
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Fig. 1. — The fraction of disks that have sufficiently low spin to 
be locally gravitationally unstable at least at one radius, as a func- 
tion of virial temperature T vlr of the halo. The computed curves 
are for z=10, but the results depend only very weakly on redshift. 
The solid and dashed curves assume gas temperatures of 8000 K (no 
H2) and 200 K (effective H2 cooling), respectively. For both tem- 
peratures, results are shown for three different disk mass-fractions, 
fd = (M disk /M halo )/(n b /f2 m ) = 0.25,0.5, and 1. For f d = 0.25, 
(as inferred for low redshift dwarfs), the majority of disks in low 
T v ; r halos will be gravitationally stable. 

^From equation (||), it is clear that if the gas cools 
from ~ 8000K to ~ 200K, this lowers by a factor of 
~ 6 the critical surface density required for local gravi- 
tational instability. In Figure pi we show the fraction of 
disk galaxies in halos with virial temperature T v ; r that 
are gravitationally unstable at redshift z = 10, for dif- 
ferent values of fd- The most important result seen in 
this figure is that in the absence of H2, i.e., at a tempera- 
ture of 8000K, virtually all halos just above the threshold 
T v ; r w 10 4 K harbor disks that are stable everywhere. In- 
stability could be promoted if a large fraction, /<j > 0.5, 
of the baryons make up the disk, however, this is unlikely 
to be the case (see discussion below). Only in relatively 
large halos (T v ; r J> 10 5 K) would most disks be able to de- 
velop instabilities. This can be easily understood with a 
simple order of magnitude calculation assuming an expo- 
nential disk which reaches the asymptotic circular veloc- 
ity on a disk scale length Rd w Ar 2 oo/"\/2, which gives: 
Q « 1.3(T gas /8000K) 1 / 2 (r vir /2 x 10 4 K)" 1 / 2 (a more de- 
tailed calculation gives a slightly lower value) . Apart from 
a weak dependence on the redshift-dependent NFW con- 
centration parameter c (at z ~ 10, we have c ~ 5 for halos 
with T V ir ~ f° w x 10 4 K), the Toomre parameter Q is in- 



dependent of redshift for a halo of given r v ; r and fd , and 
thus the results of Figure [j] depend very weakly on the 
assumed redshift. 

The figure also shows that decreasing the gas tempera- 
ture significantly increases the fraction of unstable disks, 
particularly in low mass halos. To illustrate the effect 
of the gas temperature more clearly, we have computed 
the mass-weighted fraction of disks, embedded in all halos 
with T v i r > 10 4 K, that are unstable to star formation: 



/M (Tvir = 10 *K,,) dM^(M, z)Mf uast (M, z) 
Jm(t vii =io*k,z) dM lw( M > Z ) M 

(8) 

where the mass function dN/dM(M, z) is obtained from 
standard Press-Schechter theory, and f nns t(M,z) is the 
quantity computed in Figure |l| Thus, f nas t is the quantity 
directly relevant to most cosmogonic studies, representing 
the fraction of the total gas mass reservoir that will be 
available for star-formation. The resulting unstable frac- 
tion /unst is shown as a function of redshift in Figure 0. 

The figure reveals that depending on the value of /d, 
the temperature drop caused by H2 formation and cooling 
can significantly increase (by between one or two orders of 
magnitude) the mass fraction of collapsed halos that are 
able to form stars. If H2 cooling does not proceed effi- 
ciently, then significant star formation, and therefore the 
epoch of reionization, has to likely await lower redshifts 
when still more massive (T v i r ^> 10 4 K) halos collapse, or 
until the dispersal of metals from low-spin, star-forming 
halos sufficiently contaminates surrounding halos so that 
they can cool via metal lines. 

Since £ oc A/disk, the question of whether a disk will 
be gravitationally unstable depends strongly on the mass 
fraction of baryons in the disk. As clearly revealed in Fig- 
ures [l] and ||, if fd is high, a larger fraction of disks will 
be gravitationally unstable. For simplicity, here we have 
assumed fd to be a constant, but in reality it is likely to 
vary as a function of virial temperature. One might indeed 
expect that fd is lower in halos with lower virial temper- 
atures, owing to the increased efficiency of feedback pro- 
cesses in shallow potential wells. This would accentuate 
the trend for disks in low mass halos to be gravitationally 
stable. For reference, fd ~ 0.4 yields good agreement with 
the observed sizes of disks and rotation curves over a wide 
range of circular velocities (MMW). More appropriately 
for our purposes, van den Bosch, Burkert & Swaters (2001) 
fitted angular momentum models to a sample of observed 
low-mass disk galaxies with V c i rc ~ 75 — 150 km s -1 (still 
somewhat more massive than the galaxies we are consider- 
ing) ; they found best-fit mass fractions to be considerably 
lower than the universal baryonic fraction, with a mean 
value of fd ~ 30%. Interestingly, they find that although 
disks form out of only a small fraction of the available 
baryons, they nonetheless draw most of the available an- 
gular momentum. If disks formed only out of low angular 
momentum baryons in the halo, they would be exceedingly 
compact, obviating our conclusions. 
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Fig. 2. — Mass— weighted unstable disk fraction as a function 
of redshift, as given by equation (H). At every redshift, the mass 
fraction in Toomre unstable disks have been summed over all halos 
with T v j r > 10 4 K. Thus, / uns t represents the fraction of the total 
gas mass reservoir that will be available for star— formation in all 
collapsed halos. The curves are labeled as in Figure hi If H2 cooling 
is ineffective, then the gas fraction available for star formation in 
collapsed halos is considerably lower. 

It is not clear what direction /<j should take towards 
higher redshift. While Lya cooling is likely to be more 
effective at the higher densities, typical masses and virial 
temperatures are even lower than those considered by van 
den Bosch et al. (2001), metal line cooling is absent, and 
a possibly top-heavy IMF results in more intense ionizing 
radiation (Bromm, Kudritzki & Loeb 2001) once internal 
star formation is underway and is also more likely to pro- 
duce hypernovae explosions (Heger & Woosley 2001). All 
of these effects should suppress gas cooling efficiencies and 
lower /d. If /d ~ 30%, as at low redshift, then from Figure 
||, only ^few% of the gas in disks will be available for star 
formation at redshifts z > 10. 

In summary, our results in this section underscore the 
need for H2 formation in order for most disks to be un- 
stable. In the next section, we will show that a significant 
amount of H2 is indeed expected to form in virtually all 
realistic cases; while in section §f|we consider the feedback 
processes which might suppress H2 formation and cooling. 

3. CHEMISTRY AND GAS COOLING IN T vm > 10 4 K HALOS 

In this section, we follow the coupled chemical and ther- 
mal evolution of a single fluid element when no UV flux is 
present. We find that in gas cooling from above ~ 10 4 K, 
H2 forms with a fixed abundance ih 2 ~ 10 -3 , which is 
not its equilibrium value, under a wide variety of initial 
conditions. We obtain a quantitative explanation for this 
value and the reason for the universality of this abundance. 
Readers not interested in the technical details may skip the 
rest of this section. We defer discussion of the case where 
UV flux is present to §4. 

It is not immediately obvious that H2 can form in gas 



cooling down from T v ; r > 10 4 K, since destruction of H2 by 
charge exchange and collisional dissociation is very effec- 
tive in gas at temperatures T > 5000K. There is a trough 
in the equilibrium cooling curve between ~ 8000K (when 
Lya cooling is effective) and ~ 3000K (when H2 cooling 
becomes effective), and it is possible that gas could cool 
to ~ 5000K and recombine before sufficient H2 forms to 
cool the gas to lower temperatures, hanging up in the val- 
ley between the peaks of the two cooling curves (e.g., see 
figure 12 of Barkana & Loeb 2001). 

The chemistry and cooling of primordial gas cooling 
from above 10 4 K has previously been considered by in- 
tegrating the coupled rate equations in the limit of steady 
state shock waves (Shapiro & Kang 1987, Kang & Shapiro 
1992, hereafter KS). These authors find that H2 formation 
is indeed possible and that gas can cool continuously to 
~ 200K, the limiting temperature at which H2 molecules 
comes into local thermodynamic equilibrium (LTE), and 
cooling is no longer effective. The efficacy of H2 forma- 
tion is due to the large non-equilibrium abundance of elec- 
trons in gas cooling from above T > 10 4 K; the gas cools 
and forms H2 faster than it recombines. Shapiro & Kang 
(1987) and subsequent authors have noted that H2 tends 
to form with an asymptotic abundance of ~ 10~ 3 over a 
wide range of initial conditions. We confirm this result, 
explore the process of H2 formation semi-analytically, and 
demonstrate how this asymptotic abundance can be un- 
derstood in terms of timcscale arguments. Similar results 
have been obtained recently by Susa et al. (1998). 
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Fig. 3. — The temperature (top panel), H2 and electron abun- 
dance (bottom panel) evolution as a function of time for a parcel 
of gas cooling isobarically from T = 5 X 10 4 K, with no dissociating 
UV flux. Regardless of the density of the gas, the H2 abundance 
asymptotes to a universal fraction xn 2 ~ 10 — 3 . The evolution of 
the gas is self-similar, with timescales scaling as 1/n, except at high 
gas densities n > 10 4 cm — 3 , when LTE effects become important. 
Thus, for the nT=5 X 10 6 Kcm -3 curve, the temperature evolution 
deviates from self-similarity at low temperatures (high densities). 
See text for details. 
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3.1. A universal H2 abundance 

Here, we show that over a wide range of initial condi- 
tions, Hi forms in cooling gas with a universal abundance 
xn 2 ~ 10 -3 which is not the equilibrium value. 

We begin by using a full non-equilibrium chemistry code 
to solve the rate equations for gas cooling from above 10 4 
K. The code solves the coupled set of stiff equations us- 
ing the Livermore stiff solver LSODAR, assuming primor- 
dial abundances and starting from the initial conditions 
?gas = Tvir and assuming equilibrium ionization fractions 
at Xvi r as the initial condition (for T V i r = 5 x 10 4 K, x e = 1). 
The reaction network we use in our calculations is given in 
a Table in the Appendix. The rates are identical to those 
used in Haiman, Rees & Loeb (1996), unless marked oth- 
erwise. In Figure ||, we show results from the code for gas 
cooling isobarically (the results for isochoric cooling are 
very similar). For a wide range of initial densities, the H2 
abundance asymptotes to a universal value ih 2 ~ 10~ 3 . 
We also obtain this universal H2 abundance for any start- 
ing temperature T V i r > 1.2 x 10 4 K. 

Does this universal abundance simply reflect the equi- 
librium abundance of H2 at low temperatures? The equi- 
librium abundance of Tin- is given by: 

k 9 n H n c + fci4rjH 2 n e 

77, jj = 

(fcio + fc 20 )riH + (^13 + k 2 i)n-n+ + k lg n e 
^ nB_n e -k g ^ n e -k 9 

n H fcio + «H+^i3 fcio 

The various approximations indicate the most important 
terms, which we have identified directly from the non- 
equilibrium chemistry code. The last approximation holds 
for x c < 0.03, when the H2 abundance rises rapidly and 
approaches its asymptotic value. The equilibrium abun- 
dance of H2 is then given by: 

fciomjn H - 
(fcu + fci7 + fcis)n-H+ + ki 5 nn 

&17™H+ kn 

where again the last approximation holds for x c < 0.03. 
Note that in this regime xn 2 ~ kg/kn is independent of 
density or ionization fraction and depends only on the tem- 
perature. It is interesting to note that collisional dissocia- 
tion of H2, given by H2 + H — ► 3H, is unimportant in this 
regime, contrary to previous assumptions (e.g. Omukai 
2001). Instead, the destruction of H2 is governed primar- 
ily by the charge-exchange reaction, H2 + H + — > Hj + H; 
when the gas is fully ionized, H2 formation is strongly 
suppressed. This is why recombination and H2 formation 
both initially proceed on the same timescale in Figure |[ 
We now use these expressions to check if H2 remains in 
equilibrium throughout the cooling process. In the lower 
panel of Figure |^, we plot the equilibrium abundance of H2 
(as given by the full, rather than reduced expressions for 
equations || and [lC]) , against the actual computed value of 
the H2 abundance given by the chemistry code. Initially, 
the H2 abundance follows the equilibrium value. However, 
as the gas cools to lower temperatures T < 3700K, it devi- 
ates sharply from equilibrium: the H2 abundance asymp- 
totes to the value x^ ymp ~ 10~ 3 , while the equilibrium 



value continues to rise. Thus, the universal abundance 
does not reflect the equilibrium abundance of H2 at low 
temperatures, which is significantly higher. 




2 2.5 3 3.5 4 



log(T) (K) 

Fig. 4. — The evolution of the H2 abundance as a func- 
tion of temperature, starting from T > 10 4 K. Solid lines indi- 
cate the no flux case; the dashed lines (Jii/n = 10 -3 ) and dot- 
ted lines (J21/1 = 10 _1 ) illustrate-the case when a dissociating 
UV flux is present (see section § p.2| for details; J21 has units of 
10 _21 ergs _1 cm~ 2 Hz~ 1 sr — and n has units of cm _3 L_.The bot- 
tom panel shows the equilibrium (as given by equation (id) and the 
actual (computed by the full non-equilibrium chemistry code)) H2 
abundance, while the top panel shows the timescalcs for H2 forma- 
tion tformi dissociation t<ji ss > S as cooling t coo i and electron recom- 
bination t rcc . Initially at high temperatures the H2 formation and 
destruction timescales are short compared to i sys = min(t rcc , i coo i), 
and the H2 abundance is in equilibrium, tf orm = tdiss- Once 
*diss,iform > tsys at T ~ 3700K, the H2 abundance falls out of 
equilibrium, and "freezes out" at xn 2 ~ 10 -3 . When a dissociating 
UV flux is present, tdiss asymptotes to a fixed value, instead of ris- 
ing exponentially. Eventually, tdiss becomes the shortest timescale 
in the problem. At this point, the H2 abundance falls back into 
equilibrium, and decreases rapidly. 

3.2. Why is xh 2 ~ 10" 3 ? 

Here we show that the asymptotic abundance of xh 2 ~ 
10 -3 may be quantitatively understood as a result of 
freeze-out processes. 

The H 2 will no longer evolve ( "freeze-out" ) if the H 2 for- 
mation and dissociation timescales are much longer than 
any other characteristic timescale in the system: 

tforrri) tdiss tree, tcool (H) 

where t ioim = k, x ^ {t) is the H 2 formation time, 

£diss = ~ kl7 (T)nx c is the H 2 dissociation time, t cooi 
is the cooling time, and t rrr = — ^-1— is the electron recom- 
bination time. The cooling and recombination timescales 
regulate the temperature T and the ionization fraction x e , 
which in turn affect the H2 abundance. 

We plot these timescales, as extracted from the full non- 
equilibrium calculation, in the top panel of Fig. (Q). As 
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the gas cools, the dissociation timescale tdiss grows rapidly 
and the gas falls out of equilibrium. Shortly thereafter, 
tform exceeds t rC c, t CO oi, and the H2 abundance freezes out. 

The rapid growth of tdi ss is easy to understand. The 
rates for collisional dissociation ^15 and charge exchange 
fci7 (H2 dissociation is primarily due to the latter) decrease 
exponentially as the temperature drops. On other hand, 
all other timescales t{ OIm , t coo \, t rec have only a power-law 
temperature dependence. Therefore as the gas cools the 
dissociation time rises exponentially and rapidly becomes 
much longer than all other timescales. 

Let us denote the temperature at which tdi ss > t Te c, icooi 
as Tf reeze . Since tdiss depends much more sensitively on 
temperature than t TCC ,t coo i, it exceeds them in very close 
succession. We can simply solve tdi SS = ^rec to obtain: 

k l7 {T) = k A {T) T freeze = 3700K (12) 

Up to this point, the H2 abundance is in equilibrium. The 
equilibrium value of H2 at Tf reeze is: 

xlT c « fc 9 (T frcozo )/A: 17 (r fl , czc ) a 2 x lO" 3 (13) 

which is in excellent agreement with the calculated value, 
- 1 - 2 x 10" 3 . 

For T < Tf reeze , tdiss iformj due to exponential tem- 
perature dependence of tdiss discussed above. Thus, H2 
formation is more rapid than H2 destruction and xn 2 > 
^asymp^ jj oweverj jf tfmm > t ICCl t coo i for T < Tf rooze , then 
the H2 abundance will freeze out. Since tf or m ~ ^diss ~ ^rcc 
at Tf reeze , we can write: 

jg™ = ( *** UJLV U (U) 

, I asymp I rp I ^ ' 

<-rcc \ H2 / ^ freeze / 

Since xn 2 > x^ mp , we conclude that tf orm > t rec , for 
T < Tf rcozo . Similarly, we have: 



° rm oc xh 2 T 2 x c 1 (15) 

tcool 

(where we have assumed Ah 2 oc T 4 , a good approxima- 
tion for T < 3000K). As before, xh 2 > x^ mp , so the 
minimum value of tform Adiss is given by xn 2 — x^ mp . If 
we evolve the rate equations at fixed ih 2 = x n 2 

np , the 

quantity T 2 /x c increases with time. Thus, tf orm > t coo \ 
for T < Tf ree 

zc- 

Since tdiss, ^form ^recjicoolj the conditions for freeze- 
out are satisfied. The H2 follows its equilibrium abun- 
dance until Tftoczc = 3700K, whereupon it freezes out. 
Our simple estimate from timescale arguments (|I2|), ( Jl3| ) 
for Tf reeze and a; asymp agree quantitatively with the full 
non-equilibrium calculation shown in Fig (^). This gives 
us confidence that we understand the dominant physical 
processes at work. 




j 5 1 \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
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Fig. 5. — The evolution of the gas temperature against the ion- 
ization fraction. The gas always cools along the same evolutionary 
pathway (shown in bold), regardless of whether it cools isobari- 
cally or isochorically, as long as the density stays within the lim- 
its n w 0.1 — 10 3 cm - 3 . Self-similarity is broken at low and high 
density: at low densities (dot-dashed lines) cooling off the CMB 
becomes important (the tracks here are for gas at z=20); at high 
densities n > 10 3 cm — 3 (dashed lines), the H2 coolant reaches LTE 
and cooling becomes less efficient. Also shown are the tracks when a 
dissociating UV flux is present (see section W), labelled by the value 
of J21/71. The gas follows the universal track until H2 is dissociated 
and then recombines at constant temperature. 

3.3. Self- similarity 

All of the relevant timescales in this problem scale as 
oc — . Thus, their ratio is independent of density, and the 
evolution of the gas should be 'self-similar'. There are 
three important variables controlling the evolution of the 
system: (T,x Cl n). Since the ratio of the timescales con- 
trolling (T,x c ) are independent of n, the evolution of the 
gas in the (cc e , T) plane should be independent of n (note 
that self-similarity would not be observed if x c or T were 
plotted against tn, as in e.g. Shapiro & Kang (1987), since 
in this case there is one 'hidden variable' which is missing 
from the plots). 

In Figure [|, we plot T against x c as obtained by running 
the full chemistry code and find that indeed the evolution 
of the gas is self-similar. This is subject to the following 
caveats: (i) at low density and high redshift, Compton 
cooling of the gas becomes effective; i.e., the gas cools 
faster by Compton cooling than it recombines. Compton 
cooling is important in the regime t c °^ pton < t ICC , i.e. 
when 

n < 8 x 1(T 3 ( -^—] ( cm" 3 (16) 
V8000K7 V 20 / 

Since the Compton cooling timescale is independent of 
density, the self-similar scaling is broken in the low den- 
sity regime, and the gas cooling rate is enhanced relative to 
the recombination rate (see the dot-dashed curves in Fig- 
ure ||). (ii) At high densities n > 10 4 cm~ 3 , the H 2 cooling 
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rate cx n instead of oc n 2 , as the H2 levels reach LTE (Galli 
& Palla 1998). Thus, the cooling time becomes indepen- 
dent of density, rather than t coo \ oc —, again breaking the 
self-similar scaling. In this case, the recombination rate 
is enhanced relative to the gas cooling rate (see dashed 
curves in Figure 0). However, between these two bounds 
the track in (xe, T) space followed by the cooling gas (solid 
lines in Figure o) is independent of the density and indeed 
whether the gas cools isochorically or isobarically. 

This convenient fact allows for considerable ease when 
estimating timescales, since it collapses a system of two 
variables into one variable, and allows one to 'evolve' a 
system without integrating the full rate equations. For 
the convenience of the reader we provide a fitting formula 
to this track: 

logio(T) = IQ^o+^x+a^+a^+a^ . ^ > 1Q -2 

logio(T) = b D + bix + b 2 x 2 + 6 3 x 3 + 6 4 x 4 ; x e < 10~ 2 

where x = logio(x e ) and ao = 0.683, a\ = 0.379, a 2 = 
0.522, a 3 = 0.299, a 4 = 5.67 x 10~ 2 ; and b a = 4.56, h = 
0.789, 6 2 = 9.73 x lO" 2 , 63 = 6.74 x 10~ 3 , and 64 = 
2.21 x 10~ 4 . This fit is valid to within ~ 0.5% in dex 
over the range of scales shown in Figure [| 

The cooling path followed by the gas is also indepen- 
dent of the initial temperature of the gas, as long as the 
initial ionization fraction is somewhat greater than its 
value at freeze-out, x G sa 2.4 x 10 -2 , which translates into 
T > 1.2 x 10 4 K. This is because the gas recombines at 
nearly constant temperature T ~ 9000K, and thus loses 
memory of its initial temperature and ionization fraction 
(see the plateau in temperature in Fig. [|). 

We can use the self-similar scaling between temperature 
and ionization fraction to estimate the timcscalc on which 
the H2 abundance reaches its peak value. While the tem- 
perature as a function of time is difficult to estimate, the 
ionization fraction as a function of time is straightforward. 
It is given by: 



x e (t) 



1 ~t~ ^/£rec,o 



(17) 



where x = 1 is the initial electron abundance and 

treco = l/fc4X n H = 5.6 x 10 4 (t^)" 1 yr (at 8000 K). 
At the freeze-out temperature Tf rcczc = 3700K, the ion- 
ization fraction is x c ~ 2.4 x 10~ 2 . Thus, H2 builds 
up to its peak abundance on an approximate timcscalc 



t w 2 x 10 6 ( t^L-s ) 1 

V 1 cm d / 



yr- 



For runaway collapse and fragmentation, we require 



10 3 , we have: 



^dyn 



0.22 



(t, 



,-3 



-1/2/ T \ 



(18) 



where we have assumed Ajj 2 oc T 4 , a good approxima- 
tion for T < 3000K. In T vir > 10 4 K halos, gas will al- 
ways contract via atomic cooling to high densities such 
that t coo i < tdyn- This is in contrast to T v - a < 10 4 K 
halos, where only H2 cooling is available, and only the 
high-density gas at the center of the halo can fragment. 

In summary, we find that the universal abundance 
xh 2 ~ 10~ 3 is the result of a "freeze-out" process. The H 2 
abundance follows its equilibrium value until the freeze- 
out temperature Tf rcozo 3700K, at which point the H2 



formation and destruction timescales exceed the timescales 
on which the system cools and recombines. At this temper- 
ature, the H2 abundance freezes out and stops evolving. In 
the absence of dissociating flux, H2 formation and destruc- 
tion timescales, as well as gas cooling and recombination 
timescales, all depend on collisional processes and scale as 
cx 1/n. Hence, their ratio is independent of density and 
the behavior of the gas is self-similar. 

4. FEEDBACK AND SELF-REGULATING STAR FORMATION 

In the previous section, we have seen that in the ab- 
sence of feedback, H2 always forms with an abundance of 



asymp 



10 



-3 



In this section, we consider effects that 



can, in principle, modify this conclusion. There are three 
principal feedback effects that can potentially inhibit gas 
cooling and star formation in a disk: (i) UV radiation 
and supernovae explosions heating the disk back up to 
~ 10 4 K, rendering the disk Toomre stable once again (ii) 
external dissociating UV radiation preventing H2 forma- 
tion and cooling (iii) internal dissociating UV radiation. 
In addition, we consider (iv) Lya photon trapping, which 
may cause the gas to become radiation pressure dominated 
and prevent contraction to high densities. We discuss each 
of these effects in turn. 

4.1. Gas heating and photoionization 

Our assumption of an isothermal disk is obviously an 
idealization; both thermal instability and feedback from 
star formation implies that the gas should develop a multi- 
phase structure. Nonetheless, as we have seen in §0, for 
low mass disks to remain Toomre unstable, the dominant 
thermal phase must be cold (T~ 200K). It might be a 
concern that photoionizing radiation from forming stars 
could heat the disk back up to T ~ 10 4 K, stabilizing the 
disk once again. The production rate of ionizing photons 
necessary for most of the disk to be photoionized is: 



N ion « «b / n 2 c dV 



10 56 [ ^~ 



M (T g . 



5/2 



0.05 



0.5 y Vio 4 k7 V5xio 4 k, 

5 

photons s _1 (19) 



A \ / 1 + z 



10 



For a zero-metallicity Salpeter IMF with 0.1 < M*/M < 
100 when the lifetime of the starburst is r ~ 10 6 yrs, 
the stellar ionizing photon production rate is Aj on w 

10 53 ( io^m°0 ) photons s _1 , where M*, t ot is the total mass 

in stars (Tumlinson & Shull 2000). The entire disk gas 
mass would be converted into stars before this stabilizing 
mechanism becomes important. If the stellar IMF con- 
sists wholly of supermassive stars M > 100M Q , the ioniz- 
ing photon production efficiency is an order of magnitude 
greater (Bromm et al 2001c), and this feedback mechanism 
may be marginally important. A likely much more impor- 
tant feedback mechanism could be supernova explosions, 
which can heat the gas to higher temperatures, and unbind 
it from the disk. A detailed quantitative understanding of 
supernovae feedback is still elusive (for recent progress see 
Efstathiou 2000, Mac-Low & Ferrara 1999). We merely 
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note that the effect of supernova explosions is devastating 
for low mass halos, but it becomes progressively less se- 
vere for T > 10 4 K halos. The potential wells are deeper, 
the gas contracts to the center of halos to yield higher 
binding energies, and SNe are surrounded by a denser am- 
bient medium, which can radiate away the energy of the 
explosion more efficiently. 

Another potential feedback process is production of sec- 
ondary electrons by X-ray sources. This could be an im- 
portant effect if the first light sources were mini-quasars 
(Haiman, Abel & Rees 2000), or high-redshift supernovae 
were an important source of X-rays (Oh 2001). This delays 
H2 formation and cooling, since the secondary electrons 
produced increase the rate of H2 destruction by charge 
exchange processes (Kang & Shapiro f992). As the gas 
gradually cools and recombines, the electron fraction falls 
and H2 production can proceed, unless the rate of X-ray 
heating can balance gas cooling at ~ 10 4 K to achieve sta- 
ble thermal equilibrium. 

4.2. Photodissociation by an external UV field 

Here we consider the effects of an external photodissoci- 
ating UV radiation field. We show that H2 formation and 
cooling depends primarily on a single parameter, J21/71, 
which controls the peak H2 abundance and thus the min- 
imum temperature T mm a parcel of gas can cool to. The 
value of J2i/n is significantly smaller in T V j r > 10 4 K halos: 
atomic cooling allows the gas to contract to much higher 
initial densities. Thus, the effects of an external UV radi- 
ation field are much less important, even if the gas does 
not self-shield. Throughout the rest of this paper Jzi/n is 
implicitly given in units of 10 _21 ergs _1 cmHz -1 sr~ 4 . 

4.2.1. The importance of Jixjn 

The timescale for photodissociation is given by (Draine 
& Bertoldi 1996): 



— = 2.1 x lO 4 ^ 1 /^, 



hioM y r 



(20) 



where J21 is the average flux in the LW bands in units of 
10 _21 ergs~ 1 cm~ 2 Hz~ 1 sr _1 , and the factor /shield takes 
into account the effects of H2 self-shielding for a static 
clump of gas. For now we simply set /shield = 1; we will 
consider the effects of self-shielding later, at the end of this 
section. Unlike all the other timescales previously consid- 
ered, idiss is independent of density. Photodissociation 
thus breaks the self-similarity of the no-flux case, where 
all timescales scaled as l/n. Since tj/tdiss oc Jixjn (where 
i = rec, cool, form) , the cooling history of the gas should 
depend only on the parameter J2i/n. We show this in 
detail below. 

Dependence of a^ ak on J2i/n As the gas cools and 
recombines from above 10 4 K, it initially follows the no- 
flux case, where dissociation by charge exchange is the 
shortest timescale in the problem. However, as the tem- 
perature drops, photodissociation will take over. From 

.photo .charge — cxhangc .v • . i l l 

^diss = *diss ' * ms takes place when: 



At late times (low T and x e ), LW dissociation will al- 
ways dominate H 2 destruction. Instead of rising exponen- 
tially as the temperature drops, the H2 dissociation time 
remains constant (provided self-shielding is unimportant), 
becoming the shortest timescale in the problem as t coo i, t roc 
rise. Thus, the H2 abundance rises to a peak value before 
steadily declining once t > tdiss- 

The peak abundance of H2 can be estimated as fol- 
lows. The cooling gas follows the no flux evolution along 
the universal track until the condition in equation (|2l]) is 
satisfied at some branch-out temperature T\ m . Observe 
that equation (|2|) is a function of only the temperature, 
since x c (T) is given by the fit to the universal cooling 
track. There are two possibilities: (1) Tb < Jf roozo , i.e. 



peak 



10- 



'J21 

n 



1.5 x 10~ 2 exp 



3700 K 
T 



2.4 x 10" 2 



J2i/n < 1.5 x 10~ 2 . In this case, 
The H2 abundance reaches the asymptotic freeze-out value 
and falls out of equilibrium, until t rC c,^cooi increase suffi- 
ciently that tdiss < t Tec ,t coo \. At this point, H2 falls back 
into equilibrium and its abundance steadily decreases. (2) 
Tbo > Ifreezei i-e., J-2\/n > 1.5 x 10~ 2 . In this case, pho- 
todissociation becomes important before freeze-out can oc- 
cur. H2 is always in equilibrium, and we can compute the 
peak abundance simply from the equilibrium abundance 
of H 2 at T bo , i-e. xn a (x e (T bo ),T ho ). 

These expectations are clearly borne out in runs with 
the full non-equilibrium chemistry code. In Figure (^J), 
we plot the various timescales discussed above and the 
H2 abundance as a function of temperature. We con- 
sider two different flux levels: J21/V1 = 10~ 4 (dotted lines) 
and J2i/n = 10 -3 (dashed lines), which bracket the crit- 
ical value of J21/V1 = 1.5 x 10~ 2 .For J 2 i/n = 10 _1 , the 
H2 abundance never reaches freeze out but always fol- 
lows the equilibrium value, whereas for J21/TI — 10 -3 , the 
gas reaches the freeze-out value of xn 2 ~ 10 -3 and falls 
out of equilibrium, only to fall back into equilibrium when 
tdiss = rnin(t rcc , tcooi)- The H2 abundance declines rapidly 
thereafter. In Figure ||, we see that the cooling gas follows 
the universal track until an abrupt transition takes place, 
corresponding to tdiss = nfin(i rec , t CO oi), at which point the 
gas stops cooling and recombines at constant temperature. 

Dependence of T m i n on Jzi/n The results above sug- 
gest that the minimum temperature T m { n the gas can cool 
to depends on J2i/n. We have verified this by running the 
full chemistry code for a parcel of gas at different densities 
and under different flux levels, assuming isobaric cooling 
(the result is independent of the time the gas has avail- 
able to cool in as long as tdiss < *h- If ^diss > tn, then 
photodissociation is in any case unimportant). In Figure 
H we plot the minimum temperature a parcel of gas cools 
down to, T min against J^xfn. Where the cooling of the 
gas is in the self-similar regime n ~ 10 _1 — 10 4 cm -3 , the 
curves lie on top of one another regardless of the absolute 
values of J21, n, confirming that T m i n is a function only of 
one variable J2i/n. When n > 10 4 cm~ 3 then as before 
this scaling behavior is broken: the cooling time becomes 
independent of density and the gas cools less efficiently. 
In particular T m i n is no longer a function of J2i/n but 
just J21 ( see d ashed line). We explore this regime in more 
detail in §4.3. 

We can compute T m ; n analytically. Since the gas fol- 
lows the universal track until tdiss (^21) < min(t rcc , t coo i), 



(21) 
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we can solve the equation: 



^diss — min(£ rec , £cooi) 



(22) 



to obtain the temperature T m j n at which the H2 is dissoci- 
ated and the gas stops cooling (this argument only holds 
for J 21 /n < 1.5 x 10~ 2 . For J 21 /n > 1.5 x 10" 2 then 
idiss < icooi always holds. The gas can nonetheless cool by 
a small amount, which must be calculated by solving the 
energy equation) . We plot the result as points in Figure [| 
The simple estimate idiss = ^cooi shows remarkably good 
agreement with the results of the full chemistry code. 

4.2.2. The cooled gas fraction in halos 

The significance of the previous result is as follows. In 
order for the gas to cool down to a temperature T m j n , it 
must have (^ L ) < ( „ ) - t i m particular, in the pres- 
ence of a given dissociating flux J21, the gas must have 
n > n cl -i t . In halos with T V j r < 10 4 K, H2 is the only avail- 
able coolant; if initially n < n crit the gas will never cool to 
2mi n . For this reason, H2 formation and cooling is thought 
to be inhibited in T < 10 4 K halos due to the effects of 
an external UV dissociating background (Haiman, Rees & 
Loeb 1997, Ciardi et al. 2000, Haiman, Abel & Rees 2000, 
Machacek et al. 2001). However, the effects of an exter- 
nal UV background are much less severe for T > 10 4 K 
halos, primarily because atomic cooling contracts the gas 
to sufficiently high densities that the H2 formation and 
cooling timescales are shorter than the photodissociation 
time. Initially, when the gas density is low, H2 formation 
will be suppressed, and the gas will cool by atomic cool- 
ing, contracting until n > n cr i t and H2 cooling takes over. 
From Figure [g, we see that in order to cool the gas down 
to ~ 500K, we require n > n^t ~ 10 3 J21CH1 -3 . 

^From equations (|l|) and (g) we see that almost all of 
the baryons in the disk will satisfy this criterion (as the gas 
cools it will be compressed further, increasing the gas frac- 
tion which satisfies n > n CI a). By contrast, in T vlr < 10 4 K 
halos, where there is no means of initially compressing the 
gas, only a small fraction of the gas exists at n > n cr ; t . 
We can quantify this statement by assuming that the gas 
in T V i r < 10 4 K halos is initially isothermal at the virial 
temperature of the halo. The density profile of the gas 
is then given by demanding hydrostatic equilibrium in an 
NFW halo (Makino et al. 1998). We adopt this profile to 
compute the fraction of baryonic mass of a halo above a 
given number density Mb(> n). Halos which collapse at 
similar redshift are roughly self-similar, and have a fixed 
fraction of their mass above a given density, independent 
of their mass. 

In Figure 0, we compute the fraction of gas which can 
cool down to a temperature T for halos at z = 15, as a 
function of the external UV radiation field J%\. For com- 
parison, the radiation field corresponding to n 7 ionizing 
photons per baryon in the universe is J21 ~ ■^^n 7 n;,(l + 

z) 3 x 10 21 « 10n 7 (^^) 3 (where n& is the comoving baryon 
number density, and hp is the Planck constant). Only an 
extremely small fraction of the gas can cool to low tem- 
peratures T < 500K required for a reasonably small Jeans 
mass, Mj = 10 4 (T/500 K) 3 / 2 (n/10 4 c m - 3 )- 1 / 2 M Q . 
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Analytic estimate 




/ n > 10" cm- 3 



nT=5 xl0 2 -5 xlO 5 K cm- 3 



_L 



-5 

log(J 21 /n H ) 
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min a gas par- 
for nT = 5 X 



Fig. 6. — The minimum temperature T, 
eel can cool down to, as a function of J2i/n 
10 2 , 10 3 , 10 4 , 10 5 , 10 6 , 10 7 Kcm~ 3 (solid curves). These results were 
obtained by integrating the rate equations for an isobarically cool- 
ing parcel of gas. The cooling depends only on a single variable 
• ? 2l/"(c< £ C oolAdiss); this can be seen from the fact that almost all 
the curves lie on top of one another. This scaling behaviour is only 
broken when the gas approaches high densities n > 10 4 cm — 3 (i.e., 
the nT = 5 X 10 6 ,5 X 10 7 Kcm~ 3 curves) at which point the cool- 
ing time becomes independent of density and t coo lAdiss oc J21. In 
this case, the minimum temperature the gas can cool down is not 
given by Jix/n but by J21 alone (dashed line). The points repre- 
sent the analytic estimate of T m ; n as a function of J/n by setting 
*diss = *cooli m t ne range J/n < 1.5 X 10 — 2 where this estimate is 
valid, the agreement with the full calculation is remarkably good. 

A similar quantity has been computed by Machacek et 
al. (2001) in numerical simulations of halo formation and 
gas cooling. In their Figure 3, they plot the fraction of gas 
/g as that has cooled to T < 0.5T v ; r and p > pthroshoid as a 
function of T vlI and find a strong correlation between / gas 
and T v ir- The reason for this is that the cooling function 
of H2 depends strongly on T. Since the cooling function 
Ah 2 oc T 4 is much greater at higher T, n cr i t is correspond- 
ingly lower for high T v j r halos and thus / gas is higher. We 
find reasonable agreement (to within a factor of two in the 
cold mass fraction), if we plot our results in this fashion. 
Here we choose to compute the fraction which cools to a 
fixed temperature T, since ultimately it is the final temper- 
ature which determines the Jeans mass. By our criterion, 
there is no correlation of / gas with T V i r . By comparison, 
virtually all of the disk gas in halos with T V i r > 10 4 K halos 
can cool to low temperatures T < 500K, due to the high 
densities in the disk. Overall, while levels of flux compa- 
rable to that required for the universe to be fully reion- 
ized will strongly suppress cooling in T v - a < 10 4 K halos, 
they have negligible impact on H2 formation and cooling 
in T vir > 10 4 K halos. 
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Fig. 7. — The fraction of gas in a that can cool below a given 
temperature as a function of the external UV radiation field, for ha- 
los at z=15. The radiation field corresponding to 1 ionizing photon 
per baryon in the universe is marked. In halos with T v ; r < 10 4 K 
(solid lines), only a small fraction of the gas is at sufficiently high 
density to be unaffected by external radiation fields. By contrast, in 
halos where atomic cooling operates, the gas contracts to such high 
densities that virtually all of the disk gas can form H2 and cool to 
low temperatures. 



4.2.3. Self- shielding 

We now consider the role of self-shielding. We have 
previously seen that gas in disks can form H2 and cool 
for all reasonable values of the external UV radiation 
field even if we neglect self-shielding. We therefore con- 
fine ourselves to a few general remarks. Self-shielding 
is much more important in T vu > 10 4 K halos than in 
lower mass halos: after the gas has settled in the disk, 
column densities are higher by a factor A -2 ~ 400. If H2 
forms with the asymptotic abundance of xh 2 ~ 10~ 3 , then 
from equation (0), the H2 column density is: 

10 20 exp(-r/2i? d ) (&) (, ' ^ ' ^ 1 



v0.05/ 



' l+z \ 3 / 2 
. 10 / 



5xl0 4 K y 

There are two extremes cases to be considered. If the 
gas is static v/b <C 1, then the LW flux is attenuated by 

ATh. 



a factor: f s h = min 



1, 



-0.75 



< 1 (Draine & 



Bcrtoldi 1996; this fitting formula is accurate to within a 
factor of 2 for iV H2 < 5 x 10 20 cm -2 , at which point the 
LW flux is attenuated by 5 orders of magnitude). In this 
regime, self-shielding is extremely strong and the gas is im- 
pervious to external UV radiation. However, if there are 
significant velocity gradients with v/b >> 1, then the gas 
remains optically thin to LW radiation until the damping 
wings of the LW lines overlap, at N^ 2 w 10 22 cm~ 2 (Draine 
& Bertoldi 1996, Glover & Brand 2000; in particular, see 
Fig. 2 and 3 of Glover & Brand 2000, where the fraction 
of radiation in the LW bands absorbed is computed as a 
function of H2 column density. For 7Vh 2 ~ 10 20 cm~ 2 , only 
~ 10% of the radiation is absorbed). The magnitude of 
self-shielding therefore depends on the velocity field of the 



gas, which probably lies closer to the v/b ^> 1 regime, par- 
ticularly as the disk cools by H2 cooling to temperatures 
Igas "C T v ir, and thus b -C V c i TC . We have seen that gas in 
the disk is unaffected by photodissociating UV radiation 
even if self-shielding is unimportant. On the other hand, 
the H2 optical depth affects the efficiency of H2 dissocia- 
tion by internal sources; if the optical depth is high, then 
a large fraction of the energy emitted in the LW bands by 
stars goes towards dissociating H2. We now turn to H2 
photodissociation by internal UV sources. 

4.3. Photodissociation by internal UV radiation 

H2 photodissociation by internal sources is likely to 
be the dominant source of feedback. Following Glover 
& Brand (2000), we can estimate the dissociation effi- 
ciency as follows. A single 100M© metal free star produces 
A^dis ~ 10 49 photons s _1 in the 11.15-13.6 eV range. Let 
us assume that a fraction / a b s ~ 0.1 of these photons are 
absorbed in the disk (as is appropriate if Nn 2 ~ 10 20 cm~ 2 ; 
see Figures 2 & 3 in Glover & Brand 2000) and / dis - 0.2 
of excitations lead to dissociations. Then over its main 
sequence lifetime t ms ~ 3 x 10 6 yr, the star will dissoci- 
ate M gas W /dis/ab S m H2 ^Wms/^Ha ~ 10 7 M Q of gaS of all 

its H2, assuming ih 2 ~ 10~ 3 (see also Omukai & Nishi 
(1999)). This could imply a star formation efficiency as 

low as / s tar ~ 10~ 5 . 

The efficiency with which cooling and star formation can 
still proceed depends on the amount of clumping and frag- 
mentation that can take place before internal sources of ra- 
diation turn on. Dense clumps have short free-free times 
and may be able to collapse before all the H2 in them 
is dissociated. While the H2 in the center of a cooling 
clump may initially be shielded from dissociating radia- 
tion (since ./Vh 2 3> 10 14 cm -2 for most clumps), ultimately 
if tdiss = Mh 2 /Mh 2 < tdyn, where Afn 2 is the mass of H2 
in the clump and Afu 2 is the H2 dissociation rate, then the 
clump cannot collapse. 

For a uniform density clump, the requirement 
that tdyn < ^diss translates into the require- 
ment that the clump lie at a distance D > 



10 



( iVd, 3 V /2 / n, N-7/12 f M clump \-V 6 t /.„. \ 
^lO^s- 1 ^ V10 4 cm- 3 y ^ 10 4 M o j \10- 2 J 



-1/2 



io^s- 1 ^ vio 4 cm- 3 ; ^io 4 M y V 10_2 y pc 

from the nearest dissociating source. By comparison, the 

(5^5^) (tt)" 372 ^' and 

A - ■ ~ 70pc. Even if H 2 is 



virial radius is r v 



the disk scale length is i?d ~ 7^ r vir 
completely dissociated in the disk, it may still be possible 
for dense clumps to form in the halo (where the radiation 
field from stars is weaker) , either fragmenting and forming 
stars in the halo itself or falling into the disk and accreting 
gas there. 

There is a hard limit to the temperature a parcel of 
gas exposed to dissociating radiation can cool to. This 
is because at high densities n > 10 4 cm~ 3 , collisional ex- 
citation and de-excitation dominate, the gas falls into 
LTE, and the cooling time becomes independent of den- 
sity. We can find this minimum temperature by setting 
icooi(^H 2 , T) = £diss(</2i) for the universal H2 fraction 
ih 2 ~ 10 -3 . The result is shown in Figure ^| (dashed 
curve). If the dissociating flux exceeds a critical value of 
J > Xrit ~ 10 3 , then cooling below ~ 4000K is impossi- 
ble, independent of the density of the gas. Such flux levels 
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are easily attainable from internal stellar radiation fields: 
a 100 M metal free star would generate a radiation field 

Jn « 10 3 (toS^t) (2^) 2 in the LW bands. Since 
D ~ 20pc is of order the size of typical cooling clumps, it 
implies that once a single star forms within a clump, sub- 
sequent cooling and fragmentation within the clump will 
be suppressed. 

If internal photodissociation is a very efficient source of 
feedback, then star formation initiated by H2 cooling can- 
not take place in the disk, and is likely to proceed only 
in the halo (where the sources of radiation are spread 
farther apart). If the disk has sufficiently low spin to 
be Toomre unstable at ~ 10 4 K, then unstable clumps 
could form in the disk, cooling isothermally by atomic 
cooling. As long as they can cool to extremely high 
densities (n > 10 12 cm~ 3 ), then stellar-mass size clumps 
could still form. Omukai (2001) has suggested that the 
gas cools initially by Lya and two photon emission for 
n < 10 7 cm~ 3 , and then by free-bound emission from H~ 
until n ~ 10 16 cm -2 , when the gas becomes optically thick 
to continuum radiation. This cooling route does not per- 
mit the formation of H2 by three-body reactions at high 
densities when cooling from high initial temperatures. Un- 
less sufficient H2 can form initially to cool the gas down to 
T < 2000K, collisional dissociation completely overwhelms 
the three-body H2 formation rate (see rates by Galli, Palla 
& Salpeter 1982). We have found the equilibrium values 
of H2 for T > 5000K to be negligibly small, even for ex- 
tremely high densities n ^> 10 s cm~ 3 . In principle, clumps 
cooling via atomic cooling could cool to sufficiently high 
densities so as to reach sub-solar fragmentation masses 
(Omukai 2001). This leads to a somewhat counterintu- 
itive possibility: within photodissociation regions where 
J21 > 10 3 and H2 formation and cooling is suppressed, 
the fragmentation masses could be smaller than in regions 
where H2 formation and cooling can can take efficiently. 

A possible obstacle to contraction to such high densi- 
ties is radiation pressure from Lya photon scattering; we 
examine this issue in the next subsection. If the major- 
ity of disks are Toomre stable and photodissociation by 
internal UV sources is very efficient, then efficient star 
formation may have to await not just T > 10 4 K halos, 
but halos capable of self-enriching themselves in metals, 
with sufficiently deep potential wells that metals ejected 
by supernovae rain back down on the disk in a galactic 
fountain, rather than being lost to the IGM. From local 
observations, the critical circular velocity appears to be 
v ~ 130km s^ 1 (Martin 1999). For metal-free stars radia- 
tive mass loss is likely to be inefficient (Kudritzki 2000), 
and except for stars which collapse directly to black holes 
and do not eject metals, ~half the stellar mass is con- 
verted to metals (Heger & Woosley 2001). Thus, the frac- 
tion of gas processed into stars roughly corresponds to the 
metallicity enrichment, / star ~ Z. Studies of gas cool- 
ing (Hellsten & Lin 1997, Bromm et al. 2001b) suggest 
that Z ~ 1O~ 3 Z ~ 10~ 5 is the critical metallicity at 
which metal line cooling allows the gas to cool down to 
T < 100K, implying that a prior star formation efficiency 
of /star > 10 -5 is required with H2 cooling before metal 
line cooling can take over. This seems compatible (given 
the many uncertainties) with our estimate of /star ~ 10~ 5 . 



4.4. Lya photon trapping and radiation pressure 

Collapse and fragmentation can be halted if the opac- 
ity of the contracting gas rises sufficiently that the energy 
density of trapped radiation becomes high and radiation 
pressure support becomes important. The condition for 
this to occur is: 

Lt tra .p > Eund ~ — ^— (23) 
K 

where L is the luminosity of the radiation source, ttrap is 
the timescale on which photons are trapped within the col- 
lapsing cloud, and £t>ind is the binding energy of the cloud. 
If the only source of radiation is cooling radiation and the 
system is in free fall collapse, then L sw -EbindAdyn and con- 
dition ( p3| ) translates into the condition t trap > i dy n- Rees 
& Ostriker (1977) have emphasized that if the only radia- 
tion source is cooling, then even if all radiation is trapped, 
radiation pressure cannot halt an overall collapse, since 
the resulting energy density is only half that required to 
support the cloud (P gas = 2/3J7 gas , but P rad = l/3C/ rad ; in 
addition, in our case ~ half of the energy is radiated in the 
two photon continuum, rather than in the Lya line, and 
can escape more easily). Nonetheless, they note that radi- 
ation pressure could become sufficiently important so as to 
boost the Jeans mass and halt fragmentation. We now in- 
vestigate this possibility, and also consider the case when 
an internal radiation source provides additional photons 
to provide pressure support. 

When the system is isothermally contracting by atomic 
cooling at ~ 8000K and the gas is largely neutral, the dom- 
inant opacity source is Lya scattering. The optical depth 
at line center to escaping Lya photons across a contracting 

clump of mass Mj » 3x 10 5 (g^f 2 (jQ^t^V^ 2 M 

is r = 1.2 x 10 10 ( 10 4 cm--0 \Tp~EJ' wmcn is enormous. 
However, due to the finite line width the photons diffuse 
in both space and frequency, escaping from the cloud when 
they have scattered sufficiently far from line center. This 
problem has been investigated by a number of authors 
(Adams 1975, Neufeld 1990), who find: 

ttrap = 15i C ross 10 3 < T < 10 6 
ttrap = 15i C ross(^) 1/3 r 4 1/6 T > 10 6 (24) 

where t cross — R/c is the light crossing time. In addition, 
Bonilha et al. (1979) find from Monte Carlo calculations 
that dust and velocity gradients reduce the trapping time 
by a factor / dust = (1+0 .9 d - )U .uv and / v = 1+0 2 7t? i.s respec- 
tively, where S — T^^Tdust and n « 1.5 (Ap) , where AV is 
the velocity gradient across the cloud and b is the Doppler 
parameter. Since we only consider gas of almost zero 
metallicity (Z < 1O~ 3 Z ) we ignore the effects of dust; 
in addition, since we have AV ~ b and thus / v ~ 1, veloc- 
ity gradients are unimportant. From equation ( pi| ) we find 
that at n ~ 10 4 cm -3 then £ t ra P ~ 350£ croS s ~ 8 x 10 _3 £ d yn ; 
and radiation pressure is unimportant. If the clump does 
not fragment and nr 3 — const, then ttrapAdyn °c n 8 / 18 cx 
Af ~ 8/9 and thus at n ~ 5 x 10 8 cm -3 , when Mj ~ 10 3 M , 
radiation pressure could potentially become comparable 
to the thermal gas pressure. In particular, radiation pres- 
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sure could significantly impede the contraction of a mas- 
sive clump of gas. 

In reality, the gas is likely to fragment as it cools. If 
we only consider the photon diffusion time across a Jeans 
length, then for isothermal cooling Mj oc ni? 3 oc n -1 / 2 , 
i.e. Rj oc n" 1 / 2 . This yields t trap /t dyn oc n 1 / 6 oc A/" 1 / 3 . 
In this case, the Jeans mass can fall far below a solar mass 
before radiation pressure ever becomes important. 

We now consider the case where an internal ionizing 
source, such as stars or QSOs provide additional Lya pho- 
tons from their HII regions, and consider whether collapse 
and fragmentation in the (largely neutral) gas beyond the 
surrounding regions can be halted by radiation pressure. 
Consider a cooling clump of mass Mj at T = 10 4 K. In 
order for P rad = ±£/ rad « § > p thorm = n k h T, we 

require: 

N Lya « iV ion > 4 x 10 50 f^-J photons s" 1 (25) 

where we assume that each ionizing photon results in 
the production of a Lya photon. By comparison, for 
metal free stars with M* > 100M Q , we have Ni on ~ 
10 50 (io^) photons s- 1 (Tumlinson & Shull 2000, 

Bromm et al. 2001c); therefore, only ~ 10~ 3 of the mass 
of a 1O 6 M clump need fragment into such stars for radi- 
ation pressure to become important. However, this is still 
unlikely to completely prevent smaller subclumps from col- 
lapsing. The critical star formation efficiency required to 
maintain P ra d = P gas as a function of the Jeans mass is 

Thus, when Mj ~ 300 M , / star ~ 100% of the gas must 
fragment into stars to halt collapse. We conclude that Lya 
photon pressure from stars can halt the overall contrac- 
tion of a cooling cloud, but it does not maintain the Jeans 
mass above a certain mass scale. Smaller subclumps have 
shorter trapping times and require unattainably high in- 
ternal Lya photon production rates to halt collapse. How- 
ever, Lya pressure is likely to act as a feedback mechanism 
regulating the efficiency with which such clumps collapse, 
since it prevents the ambient medium from cooling and 
condensing to high densities. 

In particular, radiation pressure from stars or AGN can 
significantly retard or halt the free fall collapse of the gas; 
if cooling gas becomes radiation pressure supported, then 
entropy fluctuations will be smoothed out and further frag- 
mentation could be very inefficient. This could well imply 
star formation efficiencies of order that given by equation 
( |2^ ) . Lya photon pressure could regulate the accretion of 
gas onto a dense clump, as well as increase the scale height 
of the disk. It is therefore likely to be an important source 
of feedback until supernova explosions become dominant. 
A full quantitative investigation of the effects of Lya pho- 
ton trapping in the early universe is beyond the scope of 
this paper (previous authors have studied the effects of ra- 
diation pressure in related contexts; see, e.g., Cox 1985; 
Bithell 1990; Haehnelt 1995). 



5. CONCLUSIONS 

In this paper, we have examined the cooling of gas in 
the first halos at high redshift with T v ; r > 10 4 K, which 
are able to cool via atomic line cooling. Previous genera- 
tions of halos with T v ; r < 10 4 K, in which the only available 
coolant is H2 molecules, are expected to have their cooling 
strongly suppressed by feedback processes (Haiman, Abel 
& Rees 2000; Omukai & Nishi 1999). Thus, the gas in 
these larger halos is likely to still be of nearly primordial 
composition. Our main findings can be summarized as 
follows: 

1. The gas should cool isothermally initially and settle 
toward a rotationally supported disk at the center of 
the dark matter halo. Unless the halo has unusually 
low angular momentum or a large fraction (/<j > 0.5) 
of the gas cools to form the disk (by contrast, present 
day dwarfs typically only have /d ~ 0.3), the major- 
ity of such disks will be gravitationally stable by the 
Toomre criterion. An additional coolant is therefore 
required to lower the gas sound speed and promote 
gravitational instability. We identify this coolant as 
H 2 . 

2. In gas which has been collisionally ionized and 
heated to a temperature of T > 10 4 K, H2 molecules 
form with a universal fraction ih 2 ~ 10~ 3 , inde- 
pendent of initial density or temperature for a wide 
range of initial conditions. This universal fraction 
is not an equilibrium value. Rather, it can be un- 
derstood as the result of a 'freeze-out'. The H2 
formation and destruction timescale become long 
compared to all other timescales at a temperature 
^freeze ~ 3700K; ih 2 ~ 10~ 3 thus reflects the equi- 
librium value of H2 at Tfreeze, the temperature at 
which the H2 abundance stops evolving. Indepen- 
dence of density can be understood from the fact 
that all collisional timescales scale as oc 1/n, and 
thus the ratio of timescales is independent of density; 
in particular, the gas follows a universal track in the 
(x e ,T) plane. Independence of the initial tempera- 
ture can be understood from the fact that the gas 
recombines isothermally at ~ 9000K, losing memory 
of the initial conditions. 

3. We have examined the feedback from UV radiation 
fields and conclude that H2 photodissociation from 
an external UV background, comparable to that re- 
quired to reionize the universe, is unlikely to be im- 
portant. This is because atomic cooling contracts 
the gas to such high initial densities that the H2 for- 
mation and cooling timescales are shorter than the 
UV photodissociation timescales. This is in contrast 
to T < 10 4 K halos where only a small fraction of 
H2 at the center of the halo is at sufficiently high 
density to promote cooling within a photodissocia- 
tion time; the fraction of gas in such halos which 
cools is much lower. We argue that radiation pres- 
sure due to Lya photon trapping could conceivably 
limit the efficiency of star formation. This is because 
only a small fraction of the gas needs to be converted 
into stars for a collapsing clump to become radiation 
pressure supported. Another possible source of pres- 
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sure support, which we have not examined, arises 
from turbulent velocity fields (Padoan 1995). 

While the fact that H 2 can form out of metal-free gas 
cooling from T > 10 4 K and cool the gas to ~ 10 2 K is a ro- 
bust result, we stress that our picture of how this enables 
gas to settle into cold, gravitationally unstable disks at the 
center of halos is considerably more uncertain. At stake is 
not only a large variation in the efficiency of star formation 
in low mass halos and thus the redshift of reionization (see 
Figures [l] and ^) , but also the mass scale at which clumps 
form, which bears upon the stellar IMF and perhaps upon 
the question of supermassive black hole formation. In par- 
ticular, the degree of fragmentation of the gas is highly un- 
certain. These important questions can only be addressed 
with some degree of confidence by high-resolution numer- 



ical simulations which are able to track the detailed gas 
hydrodynamics, chemistry and cooling, paralleling the pi- 
oneering work already done for T V j r < 10 4 K halos (Abel 
et al. 2000, Bromm et al. 2001a). The main hope of this 
paper is to stimulate further work along these lines. 
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APPENDIX: REACTION RATES AND CROSS SECTIONS 



Reaction 



Rate Coefficient (cm 3 sec 1 ) 



Reference 
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(1) H + e~ -» H+ + 2e~ 5.85 x 10- n T2exp(-157809.1/T)(l + T^y 1 [1] 

(2) He + e~ -» He+ + 2e~ 2.38 x 10~ n T^exp(— 285335. 4/T)(l + X^) _1 [2] 

(3) iJe+ + e~ ^i?e++ + 2e" 5.68 x 10~ 12 T^cxp(— 631515. 0/T)(l + T 5 *) _1 [1] 

(4) Jf+ + e- ^ H + hv 8.40x 10- 11 T-^r 3 - - 2 (l+T 6 - 7 )- 1 [1] 

(5) i?e + + e~ — > i?e + /if see expression in reference [1] 

(6) Jfe+++ e- -» #e+ + hv 3.36x 10- 10 r-^r 3 - - 2 (l+T 6 - 7 )- 1 [1] 

(7) Jf + H + —f H2 + hv see expression in reference [3] 

(8) H+ + H -y H 2 + H+ 6.40 x 1(T 10 [4] 
(9*) H + e~ H~ + hv 1.4 x i(r 18 T - 928 cxp(-T/16200) [5] 

(10) H + H 2 + e~ 1.30 x lO" 9 [6] 

(11) H+ + e~ -» 2H 1.68 x 10- 8 (T/300)-° 29 [7] 

(12) H+ + H--> H 2 + H 5.00 x 10- 6 T-5 [8] 
(13*) H~ + H+^2H 5.7 x 10- 6 T-3 +6.3 x 10- 8 - 9.2 x 10- n Ti +4.4 x 10- 13 T [5] 
(14) H 2 + e~ -» H + H~ 2.70 x 10- 8 T-i cxp(-43000/T) [9] 
(15*) H 2 + H -^3H 1.067 x 10- 10 (^) 2 012 exp[-4.463(^)- 1 ]/(l + 0.2472(^) 3 - 512 ) [10] 
(16) H 2 + H 2 — » i?2 + 2 if see expression in reference [11] 
(17*) i? 2 +H+^H+ +H 3 x 10- 10 exp(-21050/T) (T 4 < 1); 1.5 x lO" 10 cxp(-14000/T) (T 4 > 1) [5] 

(18) H 2 + e" -» 2# + e- 4.38 x 10- 10 exp(-102000/T)T a35 [12] 

(19) H- + e- ^ H + 2e- 4.00 x 10- 12 Tcxp(-8750/T) [12] 

(20) if" + H ->2H + e" 5.30 x 10- 20 T 217 cxp(-8750/T) [12] 

(21) if~ + H+^H+ + e~ see expression in reference [12] 



Cross-section (cm 2 ) 



(22) iJ + hv ^ H+ + e~ 6.3 x Kr^I^Y) 4 cxp(4 - 4(tan- 1 )e/e)/[l - exp(-27r/e)]; e = yj ' - 1 [13] 
(23*) He + /u/-»ffe++e- 0.694 x lO" 18 ^) 1 - 82 + (^) 3 ' 23 ]" 1 [14] 



(24) He+ +hv^He+++e- 1.575 x 10- 18 (^X) 4 exp(4 - 4(ta n - 1 )e/e)/[l - cxp(-27r/e)]; e = y'^ - 1[13] 

(25*) H~ + hv -> + e- 7.928 x 10 5 (^ - ^ T )^~ 3 for /iz/ > /iz/ T = 0.755 eV [15] 

(26*) i?2~ + hv —t H + see expression in reference [12] 

(27) Jf 2 + hv —* H2 + e~ see expression in reference [12] 

(28) H 2 + hv -> 2H a(v) = 3 x 10~ 18 for 11.26 cV < hv < 13.6 eV [16] 



1. Cen (1992); 2. Black (1978); 3. Rawlings ct al. (1993); 4. Karpas et al. (1979); 5. Galli and Palla (1998); 6. de Jong (1972); 7. Nakashima 
et al. (1987); 8. Dalgarno & Lcpp (1987); 9. Hirasawa (1969); 10. Dove & Mandy (1986); 11. Lcpp & Shull (1983); 12. Shapiro & Kang (1987); 
13. Osterbrock (1974); 14. Haardt & Madau 1996; 15. Abel et al. (1997); 16. Allison & Dalgarno (1970) and Dalgarno & Stephens (1970). 

The table above lists the chemical reactions included in our chemistry network. For reference, we have marked (with a 
"*") the reactions whose adopted rates are different from the compilation in Haiman, Rees & Loeb (1996). 



